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Abstract. - For an interacting system of N electrons, we study the conditions under which 
a lattice model of size L with nearest neighbor hopping t and U/r Coulomb repulsion has the 
same ground state as a continuum model. For a fixed value of N, one gets identical results when 
the inter-electron spacing to the Bohr radius ratio r s < r*. Above r*, the persistent current 
created by an enclosed flux begins to decay and r 3 ceases to be the scaling parameter. Three 
criteria giving similar r* are proposed and checked using square lattices. 



To numerically study an interacting system of N electrons, it is convenient to use a lattice 
model of size L, which can be exactly solved when N and L are small. To compare exact 
lattice results with those obtained assuming a continuum space and hence unavoidable ap- 
proximations (truncations, perturbative expansions, variational approaches, fixed node Monte 
Carlo calculations]!],...), one needs to know the role of the lattice upon the interacting system. 
Using the inter-particle spacing to the Bohr radius ratio r s , we define in this letter the value 
r* which characterizes the onset of the lattice effects for the ground state (GS). Assuming that 
N is given and varying the other parameters of the lattice model, one finds that below r* the 
lattice GS is identical to the continuum GS, r s is the scaling parameter, and the persistent 
current / created by an enclosed flux keeps its value without interaction. Above r*, the N 
electrons are localized by N lattice sites, the GSs become different, I begins to decay, and a 
lattice physics where r s is no longer the scaling parameter takes place. 

The Hamiltonian H c describing N electrons free to move on a continuum space of dimension 
d contains one body kinetic terms, two body interaction terms plus a constant term due to 
the presence of the uniform positive background which can be added in order to have charge 
neutrality. Measuring the energies in rydbergs (IRy — me 4: /2h 2 ) and the lengths in units of 
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the radius a of a sphere (circle in 2d) which encloses on the average one electron, e and m 
being the electronic charge and mass, H c reads 

1*2 1 

^ = -^E v ' + 7 E tt—ta+c™^ (!) 

Ts i=l rs l<i<j<N |r ' ^! 

which only depends on the scaling ratio r s = a/ as when N — > oo. The Bohr radius 
as = h 2 /me 2 characterizes the scale for the quantum effects. For the GS, many electrons 
are inside the quantum volume a% when r s is small and one gets a weakly coupled Fermi 
liquid. When r s is large, the volume per electron a d is large compared to a%, and one gets an 
electron solid (Wigner crystal) with almost negligible quantum effects. Hereafter, though our 
theory could be easily extended, we restrict ourselves to polarized electrons (spinlcss fcrmions) 
at d = 2. Assuming periodic boundary conditions (BCs) for a square of size D, one can ignore 
the constant term, the electronic density n s = N/ D 2 and a = 1/ ^/7rn s . 

Let us now define a square lattice model of spacing s, size L — D/s, nearest neighbor 
hopping element t = fi 2 /(2ms 2 ) and interaction strength U = e 2 /s. The lattice Hamiltonian 
Hi reads: 

The operators cj (cj) create (annihilate) a spinless fermion at the site j and (j, j') means that 
the sum is restricted to nearest neighbors, d^y is the distance between the sites j and j' in unit 
of s. When one takes periodic BCs, a convention has to be chosen for the distance d$y. One 
possible definition of djy is: 

djj' = ymin^ -j' x \,L- \j x — j' x \) 2 + mm(\j y — j' y \,L— \j y -j' y \) 2 - (3) 

Hereafter, we refer to the corresponding repulsion as the periodic singular Coulomb 

(PSC) repulsion, since it has a cusp when the intcrparticle distance dyy has one of its co- 
ordinates equal to L/2. This cusp being unphysical, we define also the periodic regularized 
Coulomb (PRC) repulsion defined from 



. L ■ 2 \3x-j' x \n . . 2 \3y-3yU 

d\\i = — sin h sm — 4 

7r V L L 

which locally coincides with the PSC repulsion, but remains analytic for all values of djy when 
s — ► 0. The PRC repulsion is essentially equivalent to the Ewald repulsion obtained from the 
periodic repetition of the considered system. In the lattice units, lRy = U 2 /At and the ratio 
r s becomes: 



The question is to know whether UL/(2tV^N) remains a scaling parameter for Hi. The 
answer is positive if r s < r*. Below r*, the lattice GS depends on the ratio UL/t for a fixed 
value of N, and coincides with the continuum GS. Let us define three criteria giving r*. 

Criterion 1: In the limit t = 0, the TV electrons are localized on N sites and form states 
| J) = . . . cj N |0) of energy Ecoul(J)- As one turns on t, one can expect that the lattice 
becomes irrelevant as each electron ceases to be localized on a single site. In analogy with the 
problem of a single particle in a disordered lattice, one can use the criterion first proposed by 
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Anderson [2] : delocalization takes place when the hopping term t between directly coupled sites 
becomes of the order of their energy spacing AE. This criterion was extended to interacting 
systems in many different contexts: onset of quantum chaos in many body spectra [3, 4, 5] 
and in the quantum computer core [6] , quasi-particle lifetime and delocalization in Fock space 
[7, 8]). In our case, the states become delocalized in the many body basis built from the states 
| J) when the matrix element (J'\Hi\J) of the one body perturbation Hi cat coupling a state 
| J) to the "first generation" of states | J') directly coupled to it by Hi exceeds their energy 
spacing AE Cou i = E Co ui(J') - E Co ui(J)- This gives t > AE Co ui- Applying this criterion to 
the GS, one obtains r* from the condition 

t^AE C oul(t = 0), (6) 

where AEcoui(t = 0) is the first level spacing of the Coulomb system on a lattice. When t 
exceeds AEc ui{t = 0), the GS is delocalized on the J-basis, and hence on the lattice, and the 
lattice GS becomes identical to the continuum GS. 

Criterion 2: Since a continuum model is invariant under translations, the motion of the 
center of mass can be decoupled from the relative motions. Thus H c can be decomposed in 
two parts, one related to the center of mass motion which is independent of the interaction, 
while the second one contains only the relative motions and hence the interaction. This has a 
very important consequence for the persistent current / driven by an enclosed Aharonov-Bohm 
flux in a continuum model: / is independent of r s and keeps its non interacting value. For 
having the topology of a 2d torus enclosing (f> along the x-direction, one takes the corresponding 
curled BC in this direction, keeping periodic BC in the y-direction. For a sufficient r s , the 
electrons form a Wigner solid and the small relative motions cannot feel the BCs. In this limit, 
/ is just given by the center of mass motion, which is independent of r s , and hence coincides 
with its non- interacting value. This point remains correct for small r s , as it was proven for 
ld-rings [9, 10, 11] and observed for d = 2 [9]. In contrast, since the previous decomposition 
into two parts does not necessary hold for Hi, I ^ I(r s = 0) for a lattice when Hi and H c 
have different GSs. The decay of / above r* (small t/U at fixed N and L) can be evaluated 
[12, 13, 14] by the leading contribution (of order N) lj; N \t/U) oc t(t/U) N ~ 1 of a t/U lattice 
expansion. The value of r s for which 

I{U = 0)*lW(t/U) (7) 

gives the criterion 2 for r* (see Fig. 2 lower left). Instead of I(<f>), one can prefer to 
use the Kohn curvature Ck = d 2 E /d(j) 2 evaluated at <p — or the GS energy change 
AtE = E ((f> = 0) — Eo((f) = 7r) when the BC is twisted in the x-direction. 

Criterion 3: When t/U — > 0, the leading correction to the Coulomb energy of Hi is ANt. 
Since the correction E v n,(r s ) to the Coulomb energy coming from the zero point vibrational 
motion of the solid when the GSs of Hi and H c are similar cannot exceed this lattice limit 
ANt (see Fig. 1), r* is also given by the condition 

E vib (r s ) w ANt, (8) 

assuming that the values of the lattice parameters can yield a Wigner solid for r s < r* . 

Let us calculate the quantities used in those criteria for the case of N — 3 spinless fermions 
on a square lattice with either the PSC or PRC repulsions. A study of the case N = 2 can 
be found in Ref.[15]. For t — 0, as shown in the insets of Figs. 1 and 2, the configurations 
minimizing the PRC and PSC Coulomb energies are different. Moving one particle by a single 
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hop increases these energies by an amount 

AF (PRC) _ 7V2ir 3 U af (psc) _ V2U 

^Coui ~ 12v /3^ 3 ' Coul 1 5 " ^ 

respectively when L is sufficiently large. For U = 0, the GS energy is given by i?o(0) = 
12t — 8t — 4icos(27r/L) for periodic BCs and becomes E (n) = 12t — 8i cos(7r/i) — At cos(3n/L) 
when one twists the BC in the x-direction. When t/U is small, A^-Eo can be calculated at 
the leading order of a i/[/-expansion [14] for N = 3. This gives when L is large: 

IA-K 2 t 9n 2 t 3 
fim A T £ « ; ^lim A T E « (10) 

where AEcoui is given by the Eqs. (9). Using these expressions, one obtains from the two 
first criteria: 

r* s (L) = AL a (11) 

where a = 4 for the PRC repulsion and a = 3 for the PSC repulsion, the constant A slightly 
depending on the taken criterion. 

For using the third criterion, one needs the zero point vibrational energy of the Wigner 
molecule that the three particles form for a sufficiently large r s < r*. Since the GSs of Hi 
and H c are identical for r s < r*(L), instead of Hi, one can use H c which is the sum of two 
decoupled terms. Denoting R = rj)/3 the coordinate of the center of mass, the first term 
reads Hqm = (^■ 2 /6'ti)V^ and corresponds to the rigid translation of the molecule while the 
other term contains the relative motions and the interaction. For a Wigner molecule, the 
second part can be simplified and expressed in terms of the normal coordinates suitable for 
describing the small vibrations around equilibrium. 

The PRC repulsion is harmonic around equilibrium, and the three particles form a diagonal 
chain as indicated in the inset of Fig. 1 when L/3 is integer. One gets four decoupled harmonic 
oscillators, two corresponding to a longitudinal mode of frequency lui — V20-B, the two others 
being a transverse mode of frequency w t — V&B, where B = (\/6e 2 7r)/(24D 3 TO). The zero 
point vibrational energy is then given for TV = 3 by: 

E V ib (r s , N = 3) = h(ui + cut) = 27r ^+^ (I) 1 rJ P (i 2 ) 

V18 V 71 / 

in rydbergs where (3 = 3/2 

When one takes the PSC repulsion, the three relative distances at equilibrium are precisely 
r = (L/2,L/2), r = (0,L/2) and r = (L/2,0) respectively when L is even. The potentials 
v(Sr) felt by the electrons around their equilibrium positions are singular and can be expanded 
as v(6r) w Cil^r^l + C^lfryl, where C\ and C*2 depend on the equilibrium positions and are 
oc e 2 /D 2 = U/L 2 . For a single particle in a ld-potcntial v(x) = C\x\, the GS energy e can 
be approximated by t/B 2 + CB where B is the GS extension and is given by de/dB = 0. 
This yields B oc (C/i) 1 / 3 and e oc (C/ 2 t/L 4 ) 1/3 . Since the 2d-potential v(6r) is separable, one 
eventually finds: 

^(r s ,iV = 3)ocr s -' 3 (13) 

in rydbergs where (3 = 4/3. As one can see, the PSC repulsion gives a higher exponent (3 when 
N = 3, which is inconsistent with the usual expansion [16] in powers of r s first proposed 
by Wigner. 
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Fig. 1. - Energy ratio Fn=3(L, U,t) as a function of r s = (U L)/(2£\/ ttN) given by the PRC repulsion 
for L = 6 (>), 9 (□), 12 (o), 15 (A), 18 (o). The dotted-dashed line gives the behavior 0.5764^71 
(harmonic vibrations of the Wigner molecule) and intersects the limiting dashed lines 12t/(4t — 
it cos(27r/L)) at the r*(L) corresponding to criterion 3. Inset: A GS configuration when t = 
and L = 24. 



From Eq. (12) and Eq. (13), since 4Nt in rydbergs is given by (AL 2 ) / (nr 2 s ) , one obtains 
from the third criterion similar results for r* as in Eq. (11), with a = 4 for the PRC repulsion 
and a = 3 for PSC repulsion. 

We now present numerical results obtained using the Lanczos algorithm, after having written 
Hi in terms of the creation (annihilation) operators of a single particle in a state of momentum 
k and having calculated the GSs in the sub-space of total momentum K — [14]. From the 
GS energy E (L, U, t) of K = 0, and for a given value of N, we define the dimensionless ratio 
F N (L,U,t) by: 

F N (L,U,t)=^ L > U > t }-^ L > U ; t = °\ (14) 
V ' ' ' E o (L,U = 0,t) V ' 

The results for the PRC repulsion are shown in Figure 1. For t = 0, the values of 
L = 6,9, 12, 15, 18 are commensurate with the period of the diagonal Wigner molecule shown 
in the inset, reducing the lattice effects. When F^—^(L,U,t) is plotted as a function of 
r s = (UL)/(2ty/nN), the different functions F/v =3 (L, U, t) scale without an observable lattice 
effect up to the r*{L) exactly given by criterion 3. Using Eq. (13) and E (L,U = 0, t) = 
12t — St — 4tcos(27r/L) one can see that the numerical results coincide with the analytical 
result Fn=3 = 0.5764 v / rJ valid for a continuum Wigner molecule. The function Fn=z{L, U, t) 
saturates to ANt/ Eq(L, U = 0,t) above r*(L), as indicated by the dashed lines. 

The corresponding results for the PSC repulsion are shown in Fig. 2 (upper left) for even 
values of L, where the GS is a triangular molecule shown in the inset when t/U — ► 0. Again 
the curves scale up to the onset r*(L) given by criterion 3. But, Fn =3 oc r^ 3 for intermediate 



6 



EUROPHYSICS LETTERS 





Fig. 2. - Upper left: Same as in Fig. 1 using the PSC repulsion, where the dotted-dashed line now 
gives the r^ 3 behavior due to the vibrations of the Wigner molecule. Inset: a GS configuration 
when t = and L = 24. Upper right: Energy ratios Fn(L, U,t) using the PRC repulsion for N = 3 
(L = 18 solid line) and N = 4 (L = 6 o, 8 □, 10 A) as a function of r s = UL/(2tVnN). Lower left: 
Dimensionless change At Eo{r s ) / At Eo(r s = 0) of the GS energy when the longitudinal BC is twisted 
for L = 6 (□), 9 (o) ,12 (A), 15 («), 18 (o) as a function of r s (PRC repulsion). Lower right: First 
energy spacing AEcoui yielded by the hop of one particle from the t = GS as a function of L. L~ 3 
behavior yielded by the PRC repulsion (A) for N = 30. Crossover towards the L~ 2 decay (dashed 
line) yielded by the PSC repulsion (□). Inset: One GS configuration yielded by the two repulsions 
(t = 0, N = 30 and L = 30). 



1 /2 

r s , as implied by Eq. (12), and not oc ry . 

In Fig. 2 (upper right), a small change of the scaling curve F N can be seen when a fourth 
electron is added, accompanied by the expected breakdown of the scaling behavior above r* . 
When N — > oo, Fjy should converge towards a thermodynamic limit depending only on r s . 
Unfortunately, a study of this convergence is out of reach of a numerical approach using exact 
diagonalization. 

In Fig. 2 (lower left), we illustrate the criterion 2, showing for TV = 3 the change AtE of 
the GS energy when a BC is twisted. Below r*, AtE does not depend on the interaction, 
while above r*, one gets the lattice limit given by Eq. (10). We have checked that the same 
conclusions can be drawn from a study of the Kohn curvature Ck{t s )- 

The PSC and PRC repulsions give rise to different r* (L), F/v =3 oc r^ 3 for the PSC repulsion, 
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differing from the conventional expansion in powers of ry . Docs this difference remain for 
larger values of N? Indeed the contribution of pairs ij having the coordinates of their spacings 
dij close to -D/2, and responsible for the r 2 J 3 behavior when dij is defined by Eq. (3), becomes 
a surface effect oc N compared to the bulk contribution oc N 2 of the remaining pairs, yielding 
^^Coui ~ AN/L 2 + BN 2 /L 3 , where A and B are constant. For a fixed L and increasing 
N, AE^^j — > BN 2 /L 3 and following criterion 1, the conventional r X J 2 expansion for Fn 
should valid for the PSC repulsion too. However, for N fixed and increasing L, the surface 
contribution dominates and the difference between the two repulsions remains in this limit. 
To illustrate this point, we have studied a Coulomb system of N = 30 electrons. The electron 
configuration minimizing the Coulomb energy is shown in the inset of Fig. 2 (lower right) for 
L = 30, i.e. one of the values of L for which the hexagonal electron lattice is commensurate 
with the underlying lattice. This configuration is given by the PRC and PSC repulsions. As 
shown in Fig. 2 (lower right) AE^c [ s [ n the crossover regime between the bulk regular 
cx L~ 3 and the singular surface cx L~ 2 behaviors, while one has the expected L~ 3 for AE^^ . 

We summarize the main results which we have checked using small lattice models and 
varying either U/t or L for N — 2,3,4. For r s < r*(L), the GSs of Hi and H c are identical, 
r s is the scaling parameter and I(r s ) = I(r s = 0). Above r*(L), r s ceases to be the scaling 
parameter and / decays. The following relations have been obtained independently of the 
definition of the Coulomb repulsion: r* cx I", E vih cx and AE Co ui oc UL^ 1 , with the 
relations a = 7 + 1 and a = 2/(2-/3) between the exponents. 

Our motivation was to investigate whether the study of a lattice model is relevant for a 
continuum model, fixing N and varying the lattice parameters. Another issue, of more direct 
physical relevance, is to study the role played by the existing lattice upon a 2d gas of correlated 
conduction electrons. We postpone the discussion of this question, where the lattice parameters 
are given, while N can be varied, to a following work. 

We thank X. Waintal for his help in the study of a Coulomb system of N = 30 electrons. 
One of us (ZAN) acknowledges the support of the European Community's Human Potential 
Programme under contract HPRN-CT-2000-00144 and the Hungarian Science Foundation 
OTKA T034832. 
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